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Abstract 

A tile Hamiltonian (TH) replaces the actual atomic interactions in a quasicrystal with effective interactions between 
and within tiles. We study Al-Co-Ni and Al-Co-Cu decagonal quasicrystals described as decorated Hexagon-Boat-Star 
(HBS) tiles using ab-initio methods. A dominant term in the TH counts the number of H, B and S tiles, favoring tilings 
of H and B only. In our model for Al-Co-Cu, chemical ordering of Cu and Co along tile edges defines tile edge arrowing. 
Unlike the edge arrowing of Penrose matching rules, however, the energetics for Al-Co-Cu do not force quasiperiodicity. 
Energetically favored structures resemble crystalline approximants to which the actual quasicrystalline compounds 
transform at low temperature. 
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Explaining thermodynamic stability is a funda- 
mental problem in the field of quasicrystals. Com- 
peting explanations range from energetic stabiliza- 
tion utilizing matching rules such as those that force 
quasiperiodicity in the Penrose tiling [Q , to entropic 
stabilization [|j focusing on the configurational en- 
tropy available in random tiling models Ex- 
perimental evidence so far has not unambiguously 
settled the matter, and the true situation is certainly 
more complex than either of the two extremes just 
described. With the advent of plausible atomistic 
quasicrystal models and advances in first-principles 
calculation methodology we hope further theoretical 
progress may be made in this area. 

Our approach reported here is based on a "tiling 
Hamiltonian" , in which a family of low energy atom- 
istic structures is placed in 1:1 correspondence with 
a family of tilings of the plane. The energetics 
of the tiling Hamiltonian is defined in a manner 
that closely approximates the ab-initio energetics 
of the atomistic structures. The energetics we de- 
rive proves remeniscient of Penrose "matching rules" 
(which force global quasiperiodicity in minimum- 
energy structures) but differs in crucial aspects. In- 
deed, we find that our tile Hamiltonian does not favor 
quasiperiodicity. Quasiperiodicity may occur at high 
temperatures as a result of random tiling configura- 



tional entropy. At low temperatures energy favors 
transformation to crystalline phases, which is indeed 
often observed experimentally (H |[ 0, |[ || . 
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Fig. 1. (A) HBS tiles and their decompositions to Penrose 
tiles, (b) Atomic decorations for AlCoCu. In (b), only 
TM and symmetry breaking Al atoms are shown. For 
AlCoNi both edge sites are decorated with Ni atoms. 

Penrose tiles are fat and thin rhombi (Fig |l|a). 
Edges are assigned single- and double- arrow deco- 
rations (as shown) which must match for common 
edges in adjacent tiles. Perfect quasicrystals obey 



these rules everywhere. The double-arrow matching 
rule causes rhombi to associate into hexagon (H) , 
boat (B) and star (S) shapes (with relative frequency 
H:B:S=a/5t:%/5:1), while the single-arrow rules force 
quasiperiodicity in the HBS tiling. It has been shown 
previously |l^ that plausible atomistic struc- 
tures of AlCoNi and AlCoCu can be described as 
HBS tilings decorated with atoms (Fig|l|b). Hence, 
we may consider the Penrose rhombus double-arrow 
rules to be satisfied by definition of our basic HBS 
tiles. 

In a tiling model of quasicrystals, the actual 
atomic interactions in the system Hamiltonian can 
be replaced with effective interactions between and 
within tiles |12 . The resulting tile Hamiltonian is 



a rearrangement of contributions to the actual total 
energy. In a simple atomic interaction picture (pair 
potentials for example) the relation between the ac- 
tual atomic interactions and the tile Hamiltonian is 
straightforward. It might be difficult to find the rela- 
tions between them for more complicated atomic in- 
teractions (many body potentials, or full ab-initio en- 
ergetics, for example) but it is theoretically possible. 
The tile Hamiltonian includes terms which depend 
only on the number of tiles, and includes other terms 
for tile interactions. The tile Hamiltonian greatly 
simplifies our understanding of the relationship be- 
tween structure and energy, and is a reasonable way 
to describe the tiling ensemble. 

Are Penrose single-arrow matching rules enforced 
by energetics of real materials? For a simple 
model 1 11 of Al7oCo9Ni2i in which both edge sites 
are occupied by Ni atoms there is no source of 
symmetry-breaking at short length scales able to de- 
fine an orientation of the tile edges. The energetics of 
structures based on HBS tiles decorated in this man- 
ner depends primarily on the numbers of H, B and 
S tiles. As seen in Fig. ||, certain phason flips con- 
vert an HS pair into a BB pair (or vice- versa). Pair- 
potential-based total energy calculations of these two 
structures |jll|, |l^ reveal that structure (a) contain- 
ing the BB pair is lower in energy than (b) contain- 
ing the HS pair by 0.2 eV. The physical origin of this 
energy difference lies in the number of 72° vertices, 
which drops by 1 in the transition HS BB. At 
a 72° vertex transition metal pairs are close neigh- 
bors, causing a reduction in the number of energet- 
ically favorable fl^ aluminum-transition metal near 
neighbor interactions. 

Hence we may express the tile Hamiltonian as 



(1) 



where Ns is the number of star tiles present, and 
the coefficient Es = 0.2 eV. To fully model decago- 
nal AlCoNi (indeed any decagonal phase) we should 
add into the Hamiltonian (eq. (|l])) terms represent- 
ing phason stacking disorder. Unfortunately, at this 




Fig. 2. Space can be tiled in many ways using HBS tiles. 
Both these approximants contain equal numbers of each 
atom type. Structures in (a) and (b) differ by the phason 
flip outhned in (b) with a dashed line, which converts two 
B tiles into and H and an S. 



time the magnitude (and even the sign) of this term 
is unknown [p^ . Additional corrections relating to 
the number of 144° vertices are small relative to the 
term shown [|ll|, Hence we focus our attention 
on the two-dimensional behavior defined by Hamil- 
tonian (cq. (0)). 

Monte Carlo simulations show that S tiles are in- 
frequent at T = lOOOK and completely absent in 
the lowest energy structures, which are random HE 
tilings with relative frequency H:B=l:r. A typical 
structure is illustrated in Fig. |^a. 

The situation for AlCoCu is more complicated 
than for AlCoNi, due to the chemical alternation of 
Co/Cu pairs on tile edges. Cockayne and Widom ||l^ 
suggested that tile edges could be assigned arrow di- 
rection based on their Co/Cu decorations (Fig.^)). 
The physical origin of Co/Cu chemical ordering rests 
on the status of Cu as a Noble Metal with completely 
filled d orbitals, unlike normal transition metals such 
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Fig. 3. Typical low temperature configurations, (a) Tile 
Hamiltonian (eq. (|])), or (eq. (||)) with Es > 2E^^. Wide 
gray bonds identify "zig-zag" structures, (b) Tile Hamil- 
tonian (eq. ^) with Es < 2Ezz- Star tiles are shaded 
gray for emphasis. 

as Co. Energetically, it turns out to be highly fa- 
vorable for Co/Cu pairs to orient such that the Co 
atoms are further removed from 72° vertices than Cu 
atoms. 

For consistency with Penrose matching rules, we 
thus define the arrow to point from Cu towards Co. 
When the HBS tiles are decorated consistently with 
the Penrose matching rules, all arrows point out- 
wards from 72° vertices, minimizing the energy asso- 
ciated with chemical ordering of Co/Cu. However 
other tilings (such as the random HB tiling illus- 
trated in Fig. ^) contain "zig-zag" structures. The 
middle of the three bonds in a zig-zag can never be 
oriented to point outwards from each of its 72° ver- 
tices, leading to a minimum energy cost for each zig- 
zag, Ezz- Hence we define our tile Hamiltonian 

H = EsNs + E,,N,, (2) 

where N^z is the number of zig-zags present, while 
the coefficients Eg = 0.2 eV and Ezz = 0.12 eV have 
been derived from full ab-initio calculations JTsl . 
Eq. (^) is actually a simplification of the full tile 
Hamiltonian O, H which captures accurately the 



energetics of the lowest energy structures. Although 
we assign the energy cost to the zig-zag shape, its ori- 
gin is the frustration of the central bond orientation, 
and not a feature of the shape itself. 

We have found a few special approximants contain- 
ing neither stars nor zig-zags. The simplest of these 
approximants (and the largest phason strain) covers 
space by translation of a single boat tile (see tiling 
Bl in Ref. ||l3|). The next larger of these approxi- 
mants (but the smaller phason strain) covers space 
with "lightbulb" objects (see Fig. ^) consisting of 
two boats and a hexagon. Other star- and zig-zag- 
free structures have been found that are basically 
superstructures of the lightbulb tiling. 




Fig. 4. zig-zag- and star-free "lightbulb tiling" structure, 
(a) Tiling has orthorhombic cell with lattice constants 
23.1 A X 31.8 A (dashed lines) containing inscribed 72° 
rhombus with 20 A edge length, (b) Model HREM struc- 
ture image ^] showing 20 A ring contrasts at vertices 
of 72° rhombus. 

For large quasicrystal approximants of low pha- 
son strain it appears impossible to simultaneously 
eliminate both stars and zig-zags. Were we to start 
with a phason strain-free random HB tiling, contain- 
ing N tiles [N/t"^ tiles of type H and Nt tiles of 
type B), a series of tile flips could segregate the tiles 
into a zig-zag-free lightbulb tiling adjoining a pure H 
tiling. Counting up the tile numbers, we see that the 
lightbulb tiling contains N / type B tiles and hence 
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7V/2t2 type H tiles. This leaves {1/t-1/2t^)N extra 
H tiles remaining to form a pure H tiling which con- 
tains 1 zig-zag per H tile. Accordingly, we conjecture 
(1/r — 1/2t'^)N is the minimum number of zig-zags 
possible in an HB tiling of N tiles at composition 
HiB,-. The number of S tiles present in an ideal HBS 
tiling of TV tiles total works out to (2/t - 1/t^)N, 
just twice the apparent minimum number of zig-zags. 
Indeed, we believe this may be the minimal allowed 
value of Ns + 2Nzz in zero phason strain tilings. If 
this were true, then the density of stars in a zero 
phason strain Penrose tiling is the minimum possi- 
ble density of stars in any zig-zag-free tiling. 

The lightbulb tiling illustrated in Fig. ^ exhibits 
a unit cell of a 72° rhombus with an edge length of 
2(cosy^ + cosff)L « 20 A where L = 6.4 A is 
the edge length of the HBS tihng for AlCoCu. 
Such a crystal structure appears when decagonal 
AI65C020CU15 is annealed at low temperatures. It 
is seen in HREM as a rhombic lattice of ring con- 
trasts identified as 20 A clusters. HREM images of 
the atomic structure associated with our lightbulb 
tiling (when decorated with atoms as in Fig. con- 
tain nearly complete ring contrasts. Fig 0b illus- 
trates a simple model high-resolution structure im- 
age obtained by superposing Gaussian functions 
at each atomic position with weight proportional to 
the atomic number (to do a better job of HREM 
modeling we should incorporate chemical and pha- 
son stacking disorder in our structure model and 
perform dynamical diffraction analysis of the elec- 
tron microscope imaging). Dark spots correspond to 
atomic columns and white to empty channels. This 
type of image should resemble HREM images from a 
thin sample near the Scherzer defocus. Fig. ^ bears 
a qualitative resemblence to the HREM patterns of 
low-temperature AI65C020CU15 in Ref. Thus it 
may be that our tiling Hamiltonian gives an indica- 
tion of the structure of the low temperature phase, 
and explains its appearance as driven by the need to 
eliminate star tiles and zig-zags. 

The precise values of the coefficients Eg and E^z 
in equations (0) and (^) can be questioned because 
they wre calculated with atoms placed at ideal sites. 
Their values will change if atomic relaxation is al- 
lowed, although we expect the general form of the 
tile Hamiltonians and the magnitudes and signs of 
the terms to be preserved. Small changes in chemical 
composition can lead to surprisingly large changes in 
the tile Hamiltonian by altering the atomic interac- 
tions specifically at those points where the unfavor- 
able star or zig-zag energies originate. Such an effect 
could explain why the low temperature structure ob- 
served for Al63C017.5Cu17.5Si2 (a 72° rhombus with 
a 51 A edge length 0) differs from that found 
in AI65C020CU15. In general, variation of the tile 
Hamiltonian parameters can lead to transitions such 




Fig. 5. Chevron tiling of 72° rhombi with 20 A edge 
lengths. The rhombic unit cell of dimensions 37.6 A x 
39.7 A matches the low temperature structure PDl of 
Al-Co-Ni § 



as that illustrated in Fig. ||. 

A chemistry dependence is also found in the case 
of AlCoNi, where small changes in composition lead 
to a wide array of different structure types Q . Tran- 
sitions as composition (or temperature) is changed 
may be related to changes in the values of terms 
in a tile Hamiltonian. For example, a change from 
Al7oCo9Ni2i to Al72CoiiNii7 results in CoAl pairs 
replacing NiNi pairs on tile edge sites at 72° ver- 
tices Consequently the energy cost of 72° ver- 
tices, and hence Es is reduced on average. However, 
the Co/Al pairs carry an edge orientation (similar 
to Co/Cu pairs) so we need to add a zig-zag energy 
into the Hamiltonian (|l|), resulting in a new Hamil- 
tonian like (^). Although the true low temperature 
phase at this composition is not certain, at a nearby 
composition of Al7iCoi4.5Nii4.5 the system indeed 
takes on one of two structures based on tilings by 72° 
rhombi with 20 A edge lengths ||]. One structure, 
known as PD2, has the unit cell of the lightbulb tiling 
(Fig. H). The other structure, known as PDl, pairs 
rhombi into "chevron" structures (Fig. |^) in which, 
again, both stars and zigzags may be avoided. Both 
of these structures have an 8 A periodicity in the 
stacking direction, so an additional term related to 
phason stacking faults may need to be included in 
our tile Hamiltonians. 

In conclusion, we show that an ensemble of low en- 
ergy quasicrystal and approximant structures may be 
modeled using very simple tiling Hamiltonians. The 
tile Hamiltonians representing Al-Co-Ni and Al-Co- 
Cu favor crystalline structures at low temperatures 
but may exhibit quasicrystals in equilibrium at high 
temperatures. The favored low energy crystal struc- 
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tures resemble the transformation products actually 
observed in these compounds at low temperatures. 
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